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Abstract 

A variety of torus models for the infrared emission of AGN has become available in literature 
over the last decade. This includes radiative transfer models using smooth or clumpy dust 
and hydrodynamic models. I will review the various types of models that are currently in use 
and point out their main similarities, differences, and limits when it comes to interpreting 
observations. Finally, current and future observational challenges for these models are 
discussed. 

1 Introduction: Why do we need torus models? 

The seminal work by Antonucci & Miller [2] showed that the prototypical type 2 active galactic 
nucleus (AGN) in NGC 1068 hosts a hidden type 1 nucleus in its center. However, due to 
significant line-of-sight obscuration of the optical and UV emission, the broad emission lines 
and "big-blue bump" emission are not directly visible to the observer. The authors found, 
however, that electrons in the polar region scatter the emission from the obscured nucleus 
and revealed the broad-emissions typical of a type 1 AGN in the polarized flux. In this sense 
the electrons in the polar region act as a mirror for the broad lines. 

The authors pictured the obscurer as a "very thick absorbing disk" without going into 
detail on its actual nature. It was Krolik & Begelman [33] who called it the "obscuring 
torus" since a toroidal geometry captures the essence of angle-dependent, geometrically- and 
optically-thick obscuration. Indeed, the idea that angle-dependent obscuration is the key 
parameter that determines if we see an AGN as type 1 or type 2 has been unequivocally 
confirmed by using jet directions in radio-loud AGN as independent inclination indicators 
[32]. 



1 Slides accompanying this review can be downloaded at http://tinyurl.com/donutsEaidcrumbs 
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It was clear early that the most likely origin of obscuration is dust. Not only does dust 
provide the obscuring properties, it also emits in the infrared (IR) and can be associated with 
the observed IR bumps in the spectral energy distributions (SEDs) of AGN [35]. One par- 
ticular feature that strongly pointed toward dust as the origin is a generic spectral turn-over 
from the big-blue bump to the IR bump at 1 ^m [13]. Considering thermal equilibrium, this 
reflects well the temperature at which dust sublimates, providing a parameter-independent 
way of explaining this generic feature (unlike non-thermal scenarios). The same considera- 
tions also point toward a size of the order of parsecs for the obscuring dusty structure [JJ. 
This brings to light a fundamental problem when investigating the torus: Although we are 
able to spectrally isolate the torus emission by observing in the infrared, the small sizes are 
beyond the spatial resolution capabilities of any single telescope. 

Finally, it has been shown that the IR power correlates well with direct luminosity 
indicators of the AGN, such as optical or UV broad-band luminosities or narrow and broad 
emission lines. More recently (e.g. see the contribution by D. Asmus) several studies showed 
that the AGN X-ray luminosity (supposedly originating from the inner disk corona) tightly 
correlates with the IR luminosity in the X-ray luminosity range of approximately Lx ~ 10 42 — 
10 erg/s, which suggests that the "torus" is a universal feature with similar characteristics 
in all Seyfert AGN [HI SUE]. 

2 The fundamentals of torus models 

This section will briefly outline the observational and physical framework for radiative trans- 
fer models of AGN tori. The purpose of these models is usually to simulate SEDs and images 
to reproduce infrared observations of AGN. These models are generally static with ad hoc 
assumptions of the dust distribution. A second family of models based on hydrodynamic sim- 
ulations aims at self-consistently predicting the mass distribution, and the radiative transfer 
is only done for specific stages of the dynamic evolution. The latter models will be briefly 
discussed in Sect. 13.41 

2.1 Observational and physical constraints 

In spite of the resolution problem, we can infer a number of constraints for the torus from 
indirect evidence. The basic framework of the torus - explaining the difference between type 
1 and type 2 AGN and the IR bump - requires that the torus is dusty, obscuring (i.e. optical 
depth Ty > 1), and geometrically thick. Further constraints on the parsec-scaled size can be 
inferred from the radiative equilibrium of dust. 

Beyond these global constraints, a lot of observational and theoretical evidence suggests 
that the dust within the torus is arranged in clumps instead of smoothly distributed. Indeed, 
the concept of a "clumpy torus" originates from the first theoretical considerations of the 
torus [3U ES] . The fact that the first models for the torus used a smooth distribution of the 
dust was owed to the lack of suitable computing power since radiative transfer calculations 
of an inhomogeneous medium are very time-expansive. 
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One of the arguments for a clumpy torus considers the velocity dispersion observed in 
the centers of galaxies [35J : If these were due to the thermal motion of (smoothly-distributed) 
dust and gas, then the corresponding temperatures in the range of 1000s K would be too large 
for dust to survive. This problem can be solved easily if the velocity dispersions reflect the 
random motion of clouds instead. It is also argued that the high HCN/CO ratio indicates that 
the gas is packed into small, dense clumps, since the emitting volume of the two molecules 
with very different critical densities must be similar to achieve such a ratio |61|, 162], Further 
evidence for a clumpy torus comes from long-term variability of the column density of type 
2 AGN [50J, although the sparse sampling provides only loose constraints on the location of 
the clouds responsible for these changes. 

From the theoretical side, a number of first-principle arguments have been discussed 
that make clumpy tori plausible. As can be shown by comparing the Jeans limit with the 
tidal forces in the gravitational field of the supermassive black hole, small compact clouds 
can be stable to both self-gravity and shear forces [351 E] • Clumpiness would also help with a 
more fundamental problem of the existence of the torus: The Eddington luminosity for dust 
around a black hole is at least 3 orders of magnitudes lower than the observed one for gas 
as a result of the relatively large ratio of cross section to mass (= opacity) [37]. As a result 
dust should be efficiently removed from the nuclei of active galaxies by radiation pressure. 
This problem can be resolved if the dust is packed into dense clouds (with strong coupling 
of gas and dust) so that the radiation pressure rather acts on the cloud as a whole than on 
individual dust grains [23]. 

2.2 Principles of radiative transfer 

The physical concept behind all radiative transfer models of AGN dust tori is thermal equi- 
librium of dust in the radiation field of the AGN. This means that the incoming radiation 
is absorbed, heats up the grain, and emits thermally according to its temperature. Dust is 
considered a "gray body" for which the efficiency of absorbing, scattering, and re-emitting a 
photon depends on the frequency of the photon. Therefore we can write the absorbed power 
of a dust grain (or cell that contains dust) as 

= / Qu A bs x F„ x A dust du. (1) 

incoming power J absorption efficiency incoming flux cross section 

The term Q U]3b b s F u notes the part of the incoming flux that is actually absorbed by the dust 
instead of scattered or transmitted. In a similar way we can write the emitted power as 

Lout = / Qv;abs x itB v (T) x A cm du. (2) 

outgoing power J absorption efficiency thermally-emitted flux emitting surface 

The detailed treatment and notion of these two equations depends on the exact setup of the 
model (e.g. single dust grain, dust cell, slab), and their solution in ID, 2D, or 3D usually 
invokes a variety of different methods, from analytic solutions to iterative and non-iterative 
Monte-Carlo (MC) simulations and raytracing techniques. 
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Figure 1: Setup of a typical model in a clumpy (right) and smooth (left) case. Both types 
of models share common characteristics, with only slight adjustments to the parameters (see 
text for details). The sketch points out the basic model properties (blue) and outlines some 
versions of the actual parametrization as used in various models (black). 



2.3 Model setup 

As mentioned previously, radiative transfer models of AGN tori use ad hoc prescriptions for 
the distribution of the dust in the model space. Despite differences in the parametrization 
used for clumpy and smooth models, and sometimes even among the same family of models, 
there are fundamental similarities that capture the basics of the torus geometry and mass 
content in any dimensionality, as illustrated in Fig. [TJ 

The model space is usually enclosed within an inner radius i?i n and an outer radius 
Rout- Whether the outer radius is defined in physical units or relative to i?; n depends on 
the choice of the actual model but does not make any difference. Since the dust sublimation 
radius (= inner radius) obeys a simple scaling law with luminosity, R m oc L 1 / 2 [4, 60J, it is 
also possible to refer to the AGN power output instead of an absolute scaling. As a cautionary 
note it is important to understand that -R out is a model parameter that does not (!) directly 
relate to an observed size. It is a mere outer boundary definition of the model space. On the 
contrary, choosing R on t too small may cause an artificial, physically unmotivated cut-off in 
the brightness distribution and severely influence the resulting SEDs and images (see Sect. 
4.1.4 in (25] for a discussion). 

The distribution of the dust in radial direction essentially defines at what radii the 
bulk of the dust mass is located relatively to all other radii and, therefore, sets which radii 
and corresponding temperatures contribute most to the overall emission (barring obscuration 
effects). A very common approach for the parametrization of the density distribution is a 
radial power law n(r) oc r~ a where a is the power-law index that defines compactness or 
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shallowness of the dust or dust-cloud distribution: for small a (0 < a < 1), the dust mass 
is distributed over a wide range of distances from the AGN, while for larger a the dust is 
concentrated toward R{ n . 

Since AGN unification requires geometrical thickness, a vertical distribution of the dust 
with scale height h/r of the order of unity or larger has to be defined. The actual parametriza- 
tion takes various shapes in literature, ranging from homogeneous vertical distributions with 
a cut-off height or scale-height to Gaussian or power-law distribution functions that cause a 
gradual fading of dust or cloud density with either height above the mid-plane or latitude. A 
Gaussian distribution has a physical motivation since it reflects the expected vertical density 
profile of an isothermal disk. 

Finally, the optical depth along any line-of-sight through the torus depends on the 
absolute density or dust mass in the model space. Common parameterizations use either of 
these values or define an optical depth value along a preferred (e.g. equatorial) line-of-site 
as a normalization. In case of clumpy models this optical depth is usually a combination of 
optical depth of an individual cloud and number of clouds along the preferred line-of-sight. 

3 An overview of AGN torus models 

Since the obscuring dust torus was established as a key component of the AGN phenomenon, 
a huge number of models have been presented to interpret IR observations. Over the years the 
models gained complexity by transitioning from smooth dust configurations, which are easy to 
simulate, to inhomogeneous (or clumpy) dust distributions. In this respect it should be noted 
again that this development is a result of increased computational power, not necessarily from 
new insight, because it was recognized from the beginning that the torus medium is probably 
clumpy. In the following, typical features and distinguishing characteristics of smooth and 
clumpy radiative transfer models, as well as hydrodynamical models, will be briefly discussed. 

3.1 Smooth models 

Early models for the torus used spherical, homogeneous dust distributions around the AGN 
or were limited to certain wavelength ranges [HET]- Major milestones were achieved when 
the dust distributions became truly 2-dimension in a sense that angle-dependent obscuration 
was incorporated into the models and the resulting anisotropic infrared emission could be 
simulated [13 E21 HH HH [191 EH [13 [391 HZ], including images that can be used to predict 
wavelength-dependent sizes [20] . Although smooth models are usually 2-dimensional by de- 
sign, 3D methods for the radiative transfer problem have been applied that made it easy to 
transition to a more inhomogeneous medium [54J. 

Despite the fact that the smooth models are probably not correctly representing the 
intrinsic structure of the torus, the relative simplicity of solving the radiative transfer problem 
allowed for adding complexity in fields that clumpy tori do usually not support. One of these 
additions is temperature/distance-dependent composition of the dust. Silicate and graphite 
grains have very different sublimation temperatures, ranging from 800 — 1000 K for silicates 
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Figure 2: Illustration of the range of parameter space for volume filling factor and optical 
depth occupied by the various clumpy torus models in literature (see references in Sect. 13, 2| ). 
Smooth models, with unity volume filling factor, are shown in blue for reference. 

and 1500 — 1900 K for graphites. In addition, at the same distance from the AGN small 
grains are hotter than larger grains, leading to a more complex temperature structure. Up to 
now, multi-species calculations have only been incorporated in smooth models [381 154] . More 
recently smooth models have also been used to test an alternative scenario and geometry 
of the dusty medium where the distribution of the dust around the AGN originates from 
hydrodynamic wind simulations [31] . 

3.2 Clumpy models 

Radiative transfer calculations of inhomogeneous media face the problem that strong density 
gradients with even stronger temperature gradients are present throughout the model space. 
This requires very fine grid sampling of the dense regions, playing against the considerable 
model volume that is essentially empty. First modeling attempts were based on the dust- 
shell approach, as used in early smooth models, by cutting "holes" into the sphere |53] , A 
more thorough way to approximate the average emission of a clumpy medium uses cloud 
source functions that are reconstructed from (2-dimensional) slab-based radiative transfer 
simulations, and the emerging emission through the inhomogeneous medium is calculated 
based on escape probabilities. More recently, advancements in computational power as well 
as the use of optimized methods allowed for direct Monte Carlo simulations (sometimes 
combined with ray-tracing techniques) of a clumpy torus, predicting both SEDs and images 
[El [22l E3 [251 EI] • These MC calculations also make it possible to incorporate and study the 
influence of inter-cloud medium, i.e. optically-thin dust that potentially occupies the void 
between the clouds [581. 
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Of course, there are lots of differences and assumptions in the treatment of the radiative 
transfer problem among the various models. To point out two distinguishing properties, Fig. [2] 
compares the optical depth of the torus/dust clouds versus the volume filling factor for the 
clumpy models. This shows that the definition of or degree of dumpiness depends on the 
specific method used. MC simulations tend to have higher volume filling factors (of the 
order of 0.1) corresponding to larger clouds for a given distance from the AGN, while the 
probabilistic approach requires very low volume filling factors, i.e. very small and dense 
clouds. Since the combination of cloud source function (size and density) and path through 
the medium (global optical depth) dominates the shape of the SEDs and the appearance of 
the torus in an image, we may expect quite some differences between the individual models. 

Another difference among the various clumpy models concerns the way the cloud density 
profile is treated. The emission-defining region of the cloud is the hot, directly AGN-heated 
surface. The bulk of the absorption and re-emission happens in this surface region leading 
to the highest temperature. When looking into literature, two very different treatments of 
this region emerge: With increasing optical depth, this region is either transitioning toward a 
"solid surface" (steep density gradient) [33j [58] or remains "fluffy" (smooth density gradient) 
[12\ |2"2"[ [55] . In the former case the hot emission of optically very thick clouds becomes a 
black-body in the near- /mid-IR, while in the latter case the hot-surface is dominated by 
optically-thin emission in the infrared independent of the total density of the cloud. 

3.3 Special purpose models 

Although not a separate family of models, a number of studies use simplified versions of a torus 
model to address special problems. Simplifications that are commonly utilized concern either 
the torus geometry (e.g. flat-disk approximation for a face-on torus) or the radiative transfer 
(e.g. using temperature/brightness distributions from full-scale models). The propagation of 
AGN variability through the dust medium has been a successful application of these kind of 
models, showing that the distribution is an important parameter that shapes the near-IR and 
mid-IR light curves [5| 1291 [26] . Another point-of-interest has been the effect of anisotropic 
accretion disk radiation on the IR emission and dust distribution at the inner rim of the torus 

3.4 Hydrodynamic models 

The advantage of radiative transfer models is their relative simplicity that makes it easy to 
simulate model grids. These grids can then be used to reproduce observations. However, they 
do not contain any physics. In order to study the dynamics of the gas and dust around the 
black hole in the region of the torus, a number of hydrodynamic models have been developed 
with the goal to reproduce the mass distribution self-consistently [631 EH E ESI E5j O 061 
[27[ I28j . This usually comes at the expense of the radiative transfer, making it difficult to 
model photometric observations for a wide range of parameters. On the other hand, since 
both densities and kinematics are predicted, the models provide a physical basis for studying 
molecular lines on scales of several to 10s of parsecs. 
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Although the details vary among the different models, it has been shown that the 
gas in the inner 5-10 pc around the AGN collapses to a thin disk [56} [65] • Outside of this 
inner dynamically cool disk the turbulence in the clumpy gas can provide for the necessary 
geometrically-thick obscuration. However, the exact vertical distribution and substructure 
does depend on the boundary conditions of the simulations. 

It has been realized recently that radiation pressure of the AGN on the circumnuclear 
medium, but also of the thermal radiation from the gas and dust itself, can have a huge effect 
on the distribution of the material in the inner torus region and the stability of dust clumps 
[36U57] . Consequently it is now attempted to work radiation pressure into the hydro dynamical 
solution |10|. [TT1 166] . One of the big problems is that in order to properly account for both 
external and internal radiation pressure, the radiative transfer through the medium has to be 
solved in some way, which adds huge overheads on the already time-consuming calculations. 

4 Caveats when using torus models 

Given the number of radiative transfer models that are available, they have become popular 
to model infrared observations of AGN. Although this is a very desirable task and potentially 
adds value to the interpretation, it is very important to keep in mind some very fundamental 
restrictions of the models: 

• Modeling means observing the model. As previously mentioned, there are different 
ways to setup a model and treat the radiative transfer that affect the SED. As such, 
if a best-fitting model is found for a set of photometric observations, these models are 
only relevant in the context of the model and do not (necessarily) tell us something 
about the real object. This fundamental problem originates from the fact that the 
infrared emission is essentially a weighted combination of gray-body emission. The 
weights are determined by the distribution of the dust and the intrinsic obscuration, 
and the same SED can be obtained by different configurations/combinations. Indeed, 
it is well possible that we will find a very good fit of the SED from a model that does 
not obey the basic requirements for an AGN torus (see Sect. [2}, in particular if the 
observations lack spatial information. 

• Different models can come to contradictory conclusions. The previous problem may also 
manifest in the answer we may get from two different torus models of the same family 
(i.e. smooth or clumpy). Given the various treatments of the radiative transfer problem 
and other "hidden" aspects (e.g. as described in the last paragraph of Sect. 13.2 j ). it is not 
surprising that the same set of parameters can lead to very different SEDs. One example 
is the appearance of the silicate feature in the different clumpy models presented in [U] 
(Fig. 15, bottom-left panel) and [25] (Fig. 7, top-right panel). Even when all other 
parameters are kept approximately the same, face-on views onto the torus results in 
an absorption feature in the former and an emission feature in the latter model (see 
the talk slides for a one-on-one comparison). Moreover, it was shown that two very 
different smooth and clumpy models can model the same set of observations equally 
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well, however never for the same set of parameters [IB]. As such, conclusions drawn 
from data fitting or parameter studies within the a given model should be considered 
model specific, even within the same family of models. 

• Model parameters are often degenerate. Finding the best fit to the model does not 
necessarily mean that the parameters are unique. In fact, since observational data 
are limited and the number of model parameters significant, it can be expected that 
similarly good-fitting solution can be found. One such combination of models is the 
opening angle, optical depth, and inclination that suffers from the fact that our view of 
an AGN is limited to one specific line-of-sight. A recent study that models photometric 
and spectroscopic IR data exemplifies the degree of parameter degeneracies that can be 
expected by analyzing 2-dimensional probability distribution functions of all pairwise 
combinations of parameters in one specific model [1] . 

5 Future challenges and possibilities 

Torus models that are currently in use face a lot of challenges from different sides. First, as 
Andy Lawrence felicitously summarized: "We have a lot of models, can we have a theory for 
them now?" Indeed, the connection of the ad hoc parameterizations of the radiative trans- 
fer models to more fundamental physics or results from (radiation- )hydrodynamic models is 
poorly explored. It is very likely that some of the configurations currently used are dynam- 
ically unstable (see contribution by M. Schartmann). Moreover, simulations and analytic 
models have shown that dust clumps exposed to the direct radiation of the AGN dissolve 
over time or are rather expelled than accreted |23} 157]. Although, admittedly, the hydrody- 
namic models are still being worked on to include all relevant physical processes, it may be 
beneficial to keep an eye on simulation results when setting up radiative transfer models, in 
particular in the realm that is not easily accessible or constrained by observations (e.g. dust 
cloud properties). 

The more direct challenges for current torus models arise from observations. With 
the advent of both infrared interferometry and high-quality photometric and spectroscopic 
surveys, our knowledge of the infrared emission of AGN has exploded. The Spitzer satellite 
provided mid- to far-IR data for a huge number of AGN of all classes. One of the features 
that emerged from these observations was a distinct near-IR bump in type 1 AGN, peaking 
about 3 — 5 /im; in fact, this near-IR bump has already been seen in earlier datasets but 
much less pronounced [13l SOI [H] . Quantitatively reproducing this bump using torus models 
has proven to be a challenge and it is generally modelled by incorporating a separate hot 
graphite-dust component. Although it has been suggested that this bump is caused by the 
fact that the dust sublimation radius is rather a sublimation region of different grain sizes 
and species, models that do account for multi-grain radiative transfer did not produce a 
significant bump |544 155j. It seems that there is great potential in improving the models and 
enhance our understanding of the geometry and composition of the hot dust by trying to 
reproduce this spectral feature. 

The biggest step forward for improving models can probably be made by taking into 
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account the spatial information from the growing number of IR-interferometrically-observed 
and -resolved AGN. For that, models have to calculate images and simulate the interferometric 
observations based on these images in Fourier space; relying on "size" comparisons may not 
be sufficient or even misleading [33] . Interferometry mainly provides information about the 
brightness distribution of the torus emission projected onto the plane of the sky. With 
the next generation of interferometric instruments (MATISSE, GRAVITY), additional such 
data will become available. Knowing the wavelength-dependent brightness distribution and 
comparing it to predictions by the model helps to eliminate degeneracies from flux-only data. 
Efforts to use radiative transfer models in this way have been undertaken only in the cases 
of NGC 1068, the Circinus galaxy, and NGC 3783 up to now [221 E3] • 
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